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ABSTRACT 

VO . 

We show that the sharp cutoff in the hard X-ray spectrum of NGC 4151, unusual 
for Seyfert 1 galaxies, can be reconciled with the average Seyfert 1 spectrum if we 
assume that the central source is completely hidden from our line of sight by the thick 
part of the accretion disk or by the broad emission line clouds. The observed X-ray 
radiation is produced by scattering of the Seyfert 1-type spectrum in the higher, cooler 
parts of the accretion disk corona, or in a wind. A sharp cutoff appears as a result 
of the Compton recoil effect. This model naturally explains a discrepancy regarding 
the inclination of the central source, inferred to be low (face-on) from observations 
of the iron Ka emission line, but inferred to be high on the basis of optical and UV 
I/-) ' observations. 
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52 ' 1- Introduction and Conclusions 

C3 ■ 

The brightest Seyfert galaxy in X-rays, NGC 4151, has a significantly different spectrum from 
the average Seyfert 1 spectrum. The average hard X-ray Seyfert 1 spectrum is well described by 
a power-law with exponential cutoff at energy, E c ~ 300 — 1000 keV, and a Compton reflection 
component. The spectrum of NGC 4151 has a much sharper decline in hard X-rays and no clear 
signature of a reflection component. While the X-ray spectra of both an average Seyfert 1 and 
NGC 4151 can be well described by models invoking thermal Comptonization of the soft radiation 
from the accretion disk in a hot corona, the corona in NGC 4151 is required to be much thicker 
(Thomson optical depth -pp ~ 2) and much cooler (T e ~ 40 — 50 keV) than the corona of Seyfert 
Is, for which r T ~ 0.2 - 0.3 and T e ~ 200 - 300 keV (Zdziarski et al. 1995, 1996). 
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In this Letter we argue that the intrinsic spectrum of NGC 4151 does not differ from that 
of Seyfert Is, if we assume that the direct component is hidden from our line of sight by the 
outer parts of the accretion disk or by the broad emission line region (BLR) close to the central 
source (Jourdain & Roques 1995) and that the observed X-ray radiation is due to scattering in the 
higher, cooler parts of the accretion disk corona or in a wind. The observed Compton-scattered 
component has much sharper cutoff than the intrinsic spectrum due to the Compton recoil 
effect. We show also that the primary X-ray spectrum of NGC 4151 is consistent with thermal 
Comptonization in active regions in the vicinity of a relatively cold accretion disk and that optical 
depths and temperatures of the hot plasma do not differ from those of Seyfert Is. Non-thermal 
models cannot be ruled out, as a strong annihilation line will be smeared out by scattering. 
The scattered component is further filtered through a complex absorber. The absorption clearly 
visible in the X-ray spectrum of NGC 4151 can be provided by the extended "atmosphere" of the 
accretion disk or BLR. 

The edge-on orientation of the NGC 4151 nucleus is strongly supported by the biconical 
geometry of the [OIII] A5007 region (Evans et al. 1993; Pedlar et al. 1993). The observed 
geometry requires the observer to be located outside the cone of UV radiation which photoionizes 
the oxygen. Recent observations of the profile of the iron Ka line, showing an extended luminous 
red wing and a sharp cutoff on the blue side, suggest an accretion disk viewed face-on (Yaqoob et 
al. 1995). This can be reconciled with the edge-on geometry deduced from the [OIII]A5007 image, 
if the central source is observed through the radiation scattered by electrons in an extended corona 
or wind. 

Finally, after correcting the UV and X-ray luminosities for dilution due to scattering, one 
also finds that NGC 4151 has luminosity ratios, Luv/Loiii and Lx/Loin, typical of Seyfert Is 
(Mulchaey et al. 1994). 



2. Scattering Model 

The scattering region is assumed to be situated along the axis of the accretion disk in the 
form of a cone (we call it the "scattering cone"). In order to avoid additional parameters, we 
assume that the radiation coming from the central source is scattered to our line of sight at a 
given angle. This angle was fixed at i = 65°, which is the best estimate obtained from optical and 
radio observations (Evans et al. 1993; Pedlar et al. 1993). 

For an arbitrary intrinsic radiation spectrum described by the photon number density, fi n ti{x), 
the Compton scattered component can be found using the Klein-Nishina formula: 

/scat(aM) oc r sc 1 + fi 2 + arxi(l - /(/) 2 /intr(Zl), (1) 

where x = hv/m c c 2 is the dimensionless photon energy, /j, = cosi, x\ = x/[l — x(l — /x)], and 
r sc is the Thomson optical depth of the scattering material. Here we assume that the electron 
temperature is much smaller than m e c 2 and that the bulk velocity in the scattering region is 
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small compared with the speed of light. Under such assumptions the scatterer can be treated as 
static and the Compton scattering cross-section for cold electrons can be applied. Then, for small 
photon energies (x <^ 0.1) the shape of the spectrum remains the same, but at higher energies, the 
cutoff appears much sharper than in the intrinsic spectrum due to the Compton recoil. 

3. Spectral Fitting 

We fit nearly simultaneous observations of NGC 4151 by ROSAT, Ginga, and OSSE in 1991 
June/July (Obs 1), and ASCA and OSSE in 1993 May (Obs 2). Description of the broad-band 
data can be found in Zdziarski et al. (1996). We used XSPEC v. 8.5 (Shafer, Haberl, & Arnaud 
1991) to fit the data. 

3.1. Exponentially Cutoff Power-Law 
with Compton Reflection 

We apply the scattering model described in § [2|, where the intrinsic spectrum is represented 
as a sum of an exponentially cut-off power-law and a Compton reflected component. This type of 
spectrum gives a good description of the spectra of Seyfert 1 galaxies (Nandra & Pounds 1994; 
Zdziarski et al. 1995). We use the Compton reflection model of Magdziarz & Zdziarski (1995). 
The parameters of the model are spectral energy index, a, cutoff energy, E c , and the amount 
of reflection, R. For an isotropic source illuminating a flat reflecting slab, R = 1 is expected. 
The iron line component is modeled by the "disk-line" model of Fabian et al. (1989). The inner 
and outer radii of the accretion disk, n and r Q , are fixed at 6 and 1000 Schwarzschild radii, re 
respectively. The line energy, Ep e , and the emissivity index, q (characterizing the line emissivity 
as oc r~ q ), are allowed to vary. The intrinsic spectrum corresponds to the inclination i = 0°. The 
scattered component is further transmitted through a complex absorber, which we approximate 
by a conventional dual absorber model with column densities N^, covering the whole source, and 
N^, covering a fraction Cp of the source (Weaver et al. 1994). In our case, by the source we 
mean the scattered component. We use the abundances from Morrison & McCammon (1983). As 
discussed in Warwick, Done, & Smith (1995) and Zdziarski et al. (1996), the data below 1 keV 
show a separate soft component which we model by a power-law with exponential cutoff at 100 
keV. An additional neutral absorber with column density, A|j oft , exceeding the Galactic value of 
2.1 • 10 20 enr 2 (Stark et al. 1992) covers the whole source. 

In spectral fitting we use only one reflection component with i = 0°, which corresponds to 
the radiation reflected from the accretion disk and scattered in the scattering cone. However, we 
expect the existence of another component arising, e.g., from the far wall of the obscuring medium 
if the latter has the form of an opaque torus (Ghisellini, Haardt, & Matt 1994; Krolik, Madau, 
& Zycki 1994). It is impossible to distinguish between these reflection components, and it is also 
quite difficult to separate reflected and scattered radiation at high energies as both have a sharp 
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cutoff due to Compton recoil. 

The Ginga data (Obs 1) provide good constraints on the amount of reflection, giving a best fit 
value of R = 0.30 ± 0.13, which is less than is typically found in Seyfert Is. In § |^ we give possible 
explanations for the weak reflection in NGC 4151. The ASCA data (Obs 2) do not constrain the 
amount of reflection, and the inferred values of R, a, and E c are strongly correlated. The data do 
not require a reflection component, but its presence cannot be ruled out. For Obs 2, we fix R at 
the best fit value for Obs 1. The best fit parameters are given in Table 1. 

The cutoff of the intrinsic spectrum (E c ~ 200 — 500 keV) appears at a much higher energy 
than in the models without scattering (E c ~ 80 — 140 keV, Zdziarski et al. 1996), and is in very 
good agreement with the cutoff of the average Seyfert 1 spectrum (E c = 560^40 keV) and the 
cutoff of the average spectrum of all Seyfert galaxies (E c = 320^^° keV, see Zdziarski et al. 1995). 
Both observed spectra of NGC 4151 are harder than the average Seyfert 1 spectrum. For Obs 1, 
the power-law index, a, of 0.7 (cf. Table 1) is inside the range of observed power-law indices of 
Seyfert Is (Nandra & Pounds 1994), while a ~ 0.5 in Obs 2 is much flatter. 



3.2. Two-Phase Disk-Corona Model 

As it is already mentioned, spectra of Seyfert galaxies can be well described by thermal 
Comptonization models (Haardt & Maraschi 1993; Stern et al. 1995). In two-phase model for an 
accretion disk-corona, soft (black body) radiation from the accretion disk gets Comptonized by 
hot thermal electron (-positron) coronal plasma. We consider two geometries of the corona: plane- 
parallel slab and localized active region in the form of a pill-box (cylinder, with height-to-radius 
ratio equal to unity). The method of Poutanen & Svensson (1996) to compute Comptonized 
spectra taking reflection from the cold disk into account is applied. The spectrum scattered by 
cold matter in the outer part of the corona is found from equation (|l|) where an intrinsic spectrum 
is the face-on spectrum of the accretion disk-corona system. The parameters of the model are 
temperature of the cold disk, Tbt>, coronal electron temperature, T e , and the optical depth of the 
corona, tt. We fix Tbb = 10 eV, because this parameter is purely constrained by the data. The 
same models for the iron line, absorber, and soft component as in § |3.1| are used. 

For the slab corona, the best fit parameters (T e ,Tr) = (210 keV, 0.24) and (218 keV, 0.30) for 
Obs 1 and 2, respectively, do not satisfy the energy balance between hot and cold phases (Haardt 
&; Maraschi 1993; Stern et al. 1995). The pill-box-corona model gives T e = 210 keV and tt = 0.54 
for Obs 1 (see Table 1). This is consistent with the situation where all energy is dissipated in an 
active region detached from the cold disk at approximately a height equal to half of its radius. 
In this condition, the X-ray source is photon-starved, and the covering factor (fraction of the 
reprocessed radiation returning to the active region) is ~ 0.4 (see Stern et al. 1995; Zdziarski et 
al. 1996). The best fit to Obs 2 gives T e = 290 keV and tt = 0.49. The energy balance requires 
dissipation of energy in an active region detached from the cold disk at a height comparable to 
its radius, implying a covering factor ~ 0.15. For pure pair corona, the compactness parameter, I 
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(for definition, see, e.g., Stern et al. 1995), is about 100 in both cases. The corresponding model 
spectra are shown in Figure [I]. 
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Table 1: Best Fit Parameters for 1991 June/July (Obs 1) and 1993 May (Obs 2) Observations of 
NGC 4151 



Note. — CPLR - exponentially cutoff power-law plus Compton 
reflection. Reflection model of Magdziarz & Zdziarski (1995) 
was used. 

PBC - pill-box-corona model by Poutanen & Svensson (1996). 
All errors are for A% 2 = 2.7. 



4. Discussion 

Observations of optical emission lines from the narrow line region in NGC 4151 indicate 
that the gas there is illuminated by an ionizing continuum stronger than the continuum observed 
from Earth by a factor 13 (Penston et al. 1990). NGC 4151 also has relatively small luminosity 
ratios Ljjy j Loin and Lx / Loin-, compared to Seyfert 1 galaxies (Mulchaey et al. 1994). This 
suggests that intrinsic UV and X-ray luminosities are underestimated by at least a factor of 
ten. Observations of the "ionizing cone" (Evans et al. 1993) clearly show the anisotropy of UV 
radiation in NGC 4151. Our models, in which the central source is completely covered by optically 
thick matter and the only radiation we see is the radiation scattered in the cone, easily explain 
these properties. 

Let us assume that the inner edge of the scattering cone is vq = 30 tq = 3 • 10 14 cm (for a 
black hole mass 5 • 10 7 -M©, Ulrich et al. 1984, Clavel et al. 1987). This is supported by the X-ray 
variability time-scale (Yaqoob & Warwick 1991). The scatterer can be produced by a wind from 
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the accretion disk. It can be highly ionized matter of normal composition or pure pair plasma. 
From observations we cannot distinguish between these alternatives. The ionization parameter, 
£ = L/r 2 n e , is of order 10 7 at r ~ 10 15 cm. Thus, highly ionized matter can scatter radiation 
out from the scattering cone without imprinting a spectral line signature on it. Assuming that 
the electron density decreases along the cone as n e = n®{r/ro)~ 2 and that r sc ~ 0.2, we can 
estimate the electron density at the inner edge of the cone, = 10 9 cm . The velocity of the 
wind would be approximately v ~ 0.1c, which is the terminal velocity of a particle if radiation 
and gravitational forces are of the same order and the particle is at rest at 10 tq from the 
center. The rate of pair production to feed the wind, n pa ir> should be about 10 48 s _1 . Assuming 
Lx = 3 • 10 44 ergs _1 and a pair yield of 1% (the fraction of luminosity converted into rest mass 
of pairs; see, e.g., Maciolek-Niedzwiecki, Zdziarski, & Coppi 1995), we get a pair production rate 
Wpair ~ 2 • 10 48 s" 1 . This shows that a pair wind is consistent with both the inferred intrinsic X-ray 
flux and the required scattering optical depth. Pairs are expected to have a Compton temperature 
of order 10 7 - 10 s K. 

One of the most intriguing properties of NGC 4151 is that different column densities of 
absorbing material are inferred from observations in UV (Ah ~ 10 19 — 10 21 cm -2 , Kriss et al. 
1992) and in X-rays (A H ~ 10 22 - 10 23 cm" 2 , Yaqoob et al. 1993). Coexistence of warm and cold 
absorbers can explain this discrepancy (Evans et al. 1993; Warwick et al. 1995). The observed 
UV and X-rays can be scattered at similar distances from the central source, consistent with the 
observed temporal correlations between the flux in these two frequency bands (Perola et al. 1986; 
Edelson et al. 1996). 



The Ginga data show that the reflection component is rather weak in NGC 4151 (see § 3.1, 
and Maisack & Yaqoob 1991). Gondek et al. (1996) found that the average Seyfert 1 also 
exhibits a deficit of reflection, having R = 0.67^g}2- If t ne underlying continuum is produced 
by Comptonization, then the spectral break appears at energies corresponding to the maximum 
of the second scattering order and is caused by anisotropy of the soft photons coming from the 
accretion disk (Stern et al. 1995). For a source compactness I ~ 100, the anisotropy break lies in 
the 2 — lOkeV energy range (see Fig. 3 in Stern et al. 1995, and dotted curves in our Fig. |TJ). 
The overall spectrum, being the sum of the Compton reflected spectrum and a broken power-law, 
is almost a perfect power-law in the 2 — 20keV band. This explains a deficit of reflection in 
Seyfert Is. If the intrinsic face-on spectrum has an amount of reflection (from the accretion disk) 
R ~ 0.6 — 0.7, then the scattered radiation is expected to have R ~ 0.3 — 0.4 due to the angular 
average. 

In the case of a torus-like geometry of obscuring matter, we can expect also some contribution 
due to reflection from the side of the torus opposite to the observer or from the outer part of the 
accretion disk. The actual fraction of the reflected radiation depends on the angular distribution 
of the intrinsic radiation, the exact geometry of the absorbing (reflecting) medium, dumpiness of 
the matter, etc. It can be about 15% for a torus geometry, assuming that optically thick material 
just covers the central source and the inclination angle is 65°, and less than 4% if the surface of 
the accretion disk is cone-like with a constant opening angle ~ 120°. This would correspond to 
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the amount of reflection R ~ 0.1 — 0.4 assuming r sc ~ 0.3 (note that in our model the intrinsic 
spectrum corresponds to i = 0°, and R is normalized to the scattered component). Since the data 
suggest R < 0.5, we conclude that either geometrical effects reduce the amount of reflection from 
the torus, the central source is anisotropic, and/or the reflecting (obscuring) matter is clumpy. 

The equivalent width of the iron line {EW ~ 200 — 300 eV) is in agreement with the 
predictions of the two-phase disk-corona models (Poutanen, Nagendra, & Svensson 1996). It is 
significantly higher than expected from the irradiation of cold matter by isotropic X-ray radiation 
having a power-law shape (George & Fabian 1991), due to the anisotropy of the Comptonized 
radiation. While the inclination of the central torus-like region is suggested to be high based on 
the Balmer line reverberation mapping (Maoz et al. 1991), recent observations of the iron line 
profile (Yaqoob et al. 1995) strongly suggest that we see the accretion disk almost face-on. This 
discrepancy has a simple explanation in the scattering model, as we see the central source via a 
"mirror" situated above the disk. 

The intrinsic spectrum of NGC 4151 can be produced by thermal Comptonization in localized 
active regions above a cold accretion disk. The compactness parameter is high, and therefore a 
significant fraction of hot plasma is expected to be in the form of electron-positron pairs. A pair 
model can also naturally explain the origin of the scattering material, as a pair wind. Despite 
the fact that a thermal model gives an excellent fit to observations, the presence of non-thermal 
processes cannot be ruled out. While non-thermal models predict a large flux at E ~ 500 keV due 
to pair annihilation, this flux is significantly diminished due to scattering by cold matter. 

An increasing fraction of the central source radiation is transmitted through the obscuring 
matter due to the Klein-Nishina decline of the Compton scattering cross-section at higher energies. 
Upper limits on the 7-ray flux at E ~ 400 keV (see Fig. [lj) provide some constraints on the 
Thomson optical thickness of the obscuring matter, r bsc- Assuming that the intrinsic flux is ten 
times larger than the observed one, we get r bsc ^ 6. 

A crucial test for the scattering model will be future observations of polarization in X-rays 
with the Spectrum- X-7 satellite (1997). The scattering model predicts a polarization of about 
20-30% perpendicular to the axis of the scattering cone. 
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grants from the Swedish Natural Science Research Council. 
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Fig. 1.— (a) Spectrum of NGC 4151 observed in June/July 1991 by ROSAT, Ginga, and OSSE. 
Solid curve represents the best fit model spectrum for pill-box-corona model (PBC, see Table 1). 
Dashed curve represents the intrinsic spectrum of the disk-corona system which is seen by the 
scattering medium (normalized to the scattered component). Dotted curve represents the underlying 
Comptonized spectrum (without reflection from the cold disk), which can be approximated by a 
broken power-law. The upper limits are 2a. (b) Same as (a), but for observations in May 1993 by 
ASCA and OSSE. 



